Unitarized pion-nucleon scattering within Heavy Baryon Chiral Perturbation Theory 
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By means of the Inverse Amplitude Method we unitarize the elastic pion-nucleon scattering ampli- 
tudes obtained from Heavy Baryon Chiral Perturbation Theory to 0(q 3 ). Within this approach we 
can enlarge their applicability range and generate the A(1232) resonance. We can find a reasonable 
description of the pion nucleon phase shifts with 0(q 2 ) parameters in agreement with the resonance 
saturation hypothesis. However, the uncertainties in the analysis of the low energy data as well as 
the large number of chiral parameters, which can have strong correlations, allow us to obtain very 
good fits with rather different sets of chiral constants. 
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Chiral Symmetry plays a fundamental role in the inter- 
actions between pions and nucleons. However, in order 
to go beyond current algebra or tree level calculations 
from simple models, one needs an effective low-energy 
field theory with all the QCD symmetries and a sys- 
tematic power counting. Heavy Baryon Chiral Pertur- 
bation Theory (HBChPT) is the best candidate to 
date for such a theory. For SU (2) chiral symmetry, the 
HBChPT degrees of freedom are the nucleons and pi- 
ons. The pions are the Nambu-Goldstone bosons of the 
spontaneous chiral symmetry breaking, whereas the nu- 
cleons are included as an isospin doublet. In the case 
of SU(3) symmetry, the pseudoscalar meson and baryon 
octets are required to describe the meson-baryon sector. 
There are also non-minimal formulations considering the 
baryon decuplet as a fundamental field || . 

HBChPT is built as an expansion in derivatives and 
meson masses, including all terms compatible with chi- 
ral symmetry. It follows the same philosophy as Chiral 
Perturbation Theory (ChPT) in the purely mesonic sec- 
tor However, there are two important differences be- 
tween them. First, the scale in the HBChPT expansion is 
not just the chiral symmetry breaking scale of pion loops 
A x = 47r/ w ~ 1.2 GeV, but also the mass of the nucle- 
ons tub — 1 GeV. Second, the nucleon four-momentum 
is of the same order as the expansion scale, no matter 
how small is the momentum transfer and even in the 
chiral limit Q. HBChPT circumvents this problem by 
treating ms as large compared to the external momenta 
and redefining the nucleon fields in terms of velocity- 
dependent eigenstates satisfying a massless Dirac equa- 
tion. The anti-baryon components can be integrated out 
||. Once this is done, it is possible to find a system- 
atic power counting in k/A x , k/rriB, <z/A x and q/rriB, 
where q stands for any meson mass or external momen- 
tum, generically denoted by 0{q). With the effective 
vertices up to a given order, one can calculate loop dia- 
grams. Each loop increases the order of the diagram so 
that any divergence can be absorbed by renormalising the 
coefficients of higher order operators. It is thus possible 



to obtain finite results order by order for any observable, 
but paying the price of introducing more chiral parame- 
ters. In particular, the full HBChPT Lagrangian up to 
0(q 3 ) has been given in w^j. It involves 5 unknown co- 
efficients to 0(q 2 ) and 23 more to 0(q 3 ), which have to 
be fitted to experiment. In this paper, we will concen- 
trate on the pion-nucleon scattering amplitude derived 
from SU(2) HBChPT to 0(q 3 ), which we will unitarize 
using the Inverse Amplitude Method (IAM). 

Despite its difficulty, several works related to ttN scat- 
tering within HBChPT have appeared in the literature 
. Here we will follow the notation of the first 0(q 3 ) 
complete result ||. Thus, we denote by and bi the 
0(q 2 ) and 0(q 3 ) parameters, respectively. The transla- 
tion to the notation of |?],|| can be found in . Only four 
0(q 2 ) and five 0(q 3 ) combinations of chiral parameters 
are relevant for 7T-N scattering. There are two sets of fit- 
ted parameters in the literature: in ||, the extrapolated 
threshold values of jl0| , together with the nuclear er-term 
and the Goldberger-Treiman discrepancy were taken as 
the experimental input, whereas in the S and P-wave 
phase shifts, somewhat away from threshold, were used as 
data for the fits. These phase shifts are given in 0] and 
are obtained as extrapolations from experimental data. 
Note that in |lj| it is suggested that the data in |Tl| l 
yield a too large a term when analyzed with HBChPT. 
As their authors pointed out, both procedures are subject 
to some caveats, either because of the many uncertain- 
ties of the data near threshold (see [0), or because the 
errors for the extrapolated phase shifts are clearly under- 
estimated (see ||). Concerning theoretical estimates, in 
H it was suggested that the a% could be understood from 
resonance exchange saturation. Further constraints from 
dispersive techniques can be found in |fl2|j . In general, 
there is a fairly good agreement between the Oj values, 
but that is not the case for 6, (see [jl2"||). 

We should stress that the predictions obtained within 
this framework are promising, although not as impressive 
H as those of ChPT for mesons. HBChPT is, of course, 
limited to low pion momentum and, with the presently 
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available calculations, certainly below ^^200 MeV M. 
The reason is, basically, that the convergence of HBChPT 
is rather slow. As a matter of fact, the contributions of 
the first three orders are frequently comparable. 

In order to improve this situation, one could go to the 
next order, thus dealing with many more parameters. 
We could also introduce more degrees of freedom, like 
the lightest resonances Jl3| , but that would also increase 
the number of parameters. In addition, some kind of 
unitarization should also be carried out in that case, to 
impose strict unitarity. Recently |14|] , remarkable results 
have been obtained by unitarizing with the N/D method 
the lowest order HBChPT at tree level, including explic- 
itly the A(1332) and iV*(1440) resonances. 

We propose an alternative approach. Encouraged by 
the HBChPT results and the success of unitarization in 
meson-meson scattering |l5| , [l6|] , we will unitarize the am- 
plitude without introducing additional fields. 



A. The IAM applied to 7T-N scattering. 

Unitarization is not foreign to effective theories. In 
fact, Pade approximants with very simple models are 
enough to describe the main features of 7r-N scattering 
Jl8|| . Although a systematic application within an ef- 
fective Lagrangian approach was called for, it was never 
carried out. 

Customarily ]T^ ], the data are presented in terms 
of partial waves of definite isospin /, orbital angu- 
lar momentum L and total angular momentum J, us- 
ing the spectroscopic notation L 2 i+i,2j+i (with L = 
S, P,... waves). Generically, the HBChPT tt-N partial 
waves, i, are obtained as a series in the momentum trans- 
fer and meson masses. Thus, basically, they are polyno- 
mials in the energy and mass variables (as well as loga- 
rithms from the loops, which provide the cuts and imagi- 
nary parts required by unitarity) . Such an expansion will 
never satisfy the 7r-N elastic unitarity condition 

Imt = q cm \t\ 2 => Imi = -q cm =>■ ~ = Rcj - i q cm , (1) 

with q cm the center of mass momentum of the incoming 
pion. But HBChPT satisfies unitarity perturbatively, i.e.: 



However, from eq.(|]) any unitarity elastic amplitude 
has exactly the following form: 



Im<i = Imt 2 = ; lmt 3 — q cm \ti\ 



(2) 



where t% stands for the 0(q k ) contribution to the ampli- 
tude. This is indeed the case in JS, but not in [0], where 
an additional redefinition of the nuclcon field allows to 
eliminate the {v ■ V) 2 /2m terms in the Lagrangian ||. 
We have performed an additional 1/m expansion of the 
results in Jp| in order to recover a pure expansion satis- 
fying eq.(||)T Thus, we count q cm and M as 0(e), so that 
each partial wave reads t ~ t\ + t 2 + t 3 + 0(e 4 ), where 
the subscript stands for the order e of each contribution. 
We have then checked that eq.(0) is verified. 
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Re{l/t)-iq c 
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for physical values of the energy, and below any inelas- 
tic threshold. The problem, of course, is how to obtain 
Re(l/t). For instance, setting Re(l/t) = |(? cm |cot<5 = 
— i + y q\ m we reobtain the familiar effective range ap- 
proximation, whereas by taking Ret ~ t\ we arrive at a 
Lippmann-Schwinger like equation Jl6| . A frequent crit- 
icism to unitarization is its apparent arbitrariness, al- 
though from eq.(||) we see that the difference between 
two unitarization methods is the way of approximating 
Re(l/t). Since we want to restrict our Lagrangian to in- 
clude just pions and nucleons preserving the SU(2) chiral 
symmetry, the most general approximation to Re(l/t) is 
HBChPT. In particular, we will take the 0(q 3 ) calcula- 
tions, but the method can be easily generalized to higher 
orders. Thus, we arrive at 



h-h + q/h-Reta 



(4) 



where we have only kept the relevant order in Re[(ii +t 2 + 
t3) -1 ]. This is the 0(q 3 ) form of the IAM. Note that if 
we reexpand in powers of q, we recover at low energies the 
HBChPT result. However, as it is written, the amplitude 
explicitly satisfies elastic unitarity. Furthermore, using 
eq.(H), we can rewrite Re £3 + iq cm t\ — t 3 , which can 
be analytically continued to the complex plane (where, 
for instance, we will look for the pole associated to the 
A(1232)). Incidentally, eq.([|) thus rewritten is a Pade 
approximant of the 0(q 3 ) series. From the K matrix 
point of view, we have identified K^ 1 = Ret -1 . 

Even though the elastic unitarity condition is only sat- 
isfied for real values of s above threshold, the use of the 
IAM in the complex plane can be justified using disper- 
sion relations ||15| , provided one is not very far from the 
physical cut. In other regions (around the left cut for 
instance) the IAM would be inappropriate. As a conse- 
quence, it is also possible to reproduce the poles in the 
second Riemann sheet, which are close to the physical 
cut and are associated to resonances. As a matter of fact 
the IAM has been successfully applied to meson-meson 
scattering |l5). In particular, using the 0(p 4 ) ChPT La- 
grangian, the IAM generalized to coupled channels yields 
a remarkable description of all channels up to 1.2 GeV, 
including seven resonances [ |l6| ]. Using the Lippmann- 
Schwinger like equation mentioned above it is also possi- 
ble to describe the S-wave kaon-nucleon scattering, using 
the lowest order Lagrangian fL9|] , including the A(1405). 

Let us then use eq.(|j) when t are the L2/+1, 2,7+1 
N partial waves. The resulting amplitudes will be fitted 
to the phase shifts, which are actually an extrapo- 
lation, not including the experimental errors. For the fit 
we have used the MINUIT Function Minimization and 
Error Analysis routine from the CERN program Library. 
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As it is customarily done in the literature, we will assign 
an error to the data in For instance, in ref. JtJ the 
central points have been given a 3% uncertainty. How- 
ever, since our fits will cover wide energy ranges, the use 
of a constant relative error will give more weight to the 
low energy data. Thus we have also added an additional 
systematic error of 1 degree. ( A 5% error plus a y/2 sys- 
tematic error was used in |14| .) This error is needed to 
use the minimization routine, and, although the order of 
magnitude may seem appropriate, the values are rather 
arbitrary, so that the meaning of the % 2 /d.o.f. obtained 
from MINUIT has to be interpreted cautiously. 

Furthermore, the data near threshold are subject to 
many uncertainties, so that, also following Q|, we will 
start our fits at y/s =1130 MeV. Hence the threshold 
parameters are real predictions in our approach. In ad- 
dition, we should limit the approach to energies where in- 
elasticities can be neglected. In particular, we will not use 
our Pii, P 13 and P31 phase shifts above the irwN thresh- 
old (^/s ~ 1220 MeV), since they are very small and 
inelasticities could be significant. The S31 and Su phase 
shifts are larger and we fitted them up to ~ 1360 MeV. 
The inelasticities in P33 are negligible up to 1400 MeV 
since this channel is dominated by the A(1232) which is 
strongly coupled to 7T-N. 

a. The IAM and Resonance Saturation. Following 
the suggestion that the 0(p 2 ) parameters can be under- 
stood from resonance saturation [gj , it is natural to try to 
make an IAM fit constrained with this hypothesis. Thus, 
we first fix the to the values of || , which are compatible 
with the saturation hypothesis. The resulting 0(p 3 ) pa- 
rameters are given in Table |. In addition, we give in Ta- 
ble [l| the values for a second fit where we have allowed the 
cii parameters to vary within the ranges expected from 
resonance saturation. The results of fit 2 are plotted, as 
a solid line, in Figure |[ We used "hatted" quantities, 
bi, because their values do not necessarily correspond to 
those of HBChPT since now they are also absorbing the 
IAM resummation effects and some high energy informa- 
tion. Only if there was a very good convergence of the 
theory at low energies the values of bi should be similar 
to the bi (as it happens in ChPT). Therefore, at present, 
our bi should not be used to calculate any other process 
at low-energies. 

Not surprisingly, there are strong correlations between 
parameters. Unfortunately, from MINUIT we cannot get 
the actual form of the correlation. However, looking only 
at linear combinations with integer coefficients, some of 
them, like b\ + bi + 63 or b\ + 62 + 263 + b\§ — 615, remain 
within natural sizes for these fits. Nevertheless, as we 
have commented, there is a considerable uncertainty in 
the precise values of the bi at present (see Table 2 in [|12[ ) . 

In summary, the main conclusion from Figure [l] is that 
it is possible to obtain an improved description of ttN 
scattering including the A(1232), with the ai values ob- 
tained from resonance saturation. Note however, that the 
S31 phase shift does not have any real improvement. For 



illustration we also give in Figure |l| the extrapolation of 
the 0(q 3 ) HBChPT results to high energies (dotted line) 
as well as the IAM result (dashed line), using the and 
bi values in ||. 

b. Unconstrained IAM fits. Of course, we can get 
much better fits (all with x 2 / d.o.f .<\) by leaving all the 
parameters free. For illustration, see the dashed-dotted 
line in Figure [l]. There are again strong correlations and 
the actual value of each one of the di and bi could be ex- 
tremely unnatural. The correlations now are even more 
complicated due to the quadratic t\ term in the denom- 
inator of eq.(|]). By inspection of the analytic formulas, 
we find that the a\ + &2, 65 — 4&i, 1{b\ + 62) + (big — 615) 
combinations are the most relevant, and remain rather 
stable for these fits. However it is not possible to obtain 
a meaningful determination of each individual parameter 
without any other additional assumption (like resonance 
saturation). That is again due to the large number of 
parameters, but also to the slow HBChPT convergence. 

c. The A(1232) resonance. The IAM generates di- 
namically a pole in the second Riemann sheet at yfs — 
(1212 — i 47) MeV, which is rather stable within all the 
fits and in very good agreement with the data pcfl . 

d. Threshold parameters. For definitions and nota- 
tion we refer again to [Bj . Our results are shown in Table 
||, where we have also listed the experimental values, 
extracted from As pointed out in [pi and H, the 
errors for those values are clearly underestimated. Hence, 
it should be borne in mind that the threshold parameters 
are not so well determined as it may seem from those er- 
rors (see also ) . Our fits give a reasonably good agree- 
ment with experiment for the S-wave scattering lengths. 
For most P-waves, we agree with the order of magnitude 
and sign. Our results are also in rough agreement with 
§, where they give: -0.07 GeV" 1 < a+ < 0.04 GeV" 1 
and 0.6 GeV" 1 < a^ < 0.67 GeV" 1 . 

B. Conclusions and discussion 

We have unitarized the HBChPT 0(q 3 ) ttN elas- 
tic scattering amplitude with the Inverse Amplitude 
Method. This approach is able to describe the phase 
shifts up to the inelastic thresholds and, in addition, it 
gives the correct pole for the A(1232) in the P33 chan- 
nel. Our fits use the extrapolated phase shifts between 
■s/s = 1130 MeV and the corresponding inelastic thresh- 
olds. Within this approach, we can predict the thresholds 
values, which for the S waves are in good agreement with 
experiment and with recent determinations. 

Unfortunately, since there are large correlations be- 
tween some parameters, it is possible to obtain good fits 
with very different sets of parameters, which can have 
rather unphysical values. This is due to different reasons: 
a) The slow convergence of the series, since contributions 
from different orders are comparable in almost every par- 
tial wave. The effect of higher order terms, which was 
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less relevant at threshold, is absorbed in our case in the 
values of the chiral coefficients, b) There are strong cor- 
relations between the parameters, and the fits are only 
sensitive to certain combinations. Hence the values of 
each individual coefficient are meaningless. 

The most relevant conclusion of this study is that we 
can still reproduce the A(1232) with the at values ex- 
pected from the resonance saturation hypothesis, keeping 
a reasonably good description for the other channels. 

Finally, we would like to remark that the method de- 
veloped here can be easily extended to the case of SU (3) 
symmetry as well as to the coupled channel formalism. 
Further work along these lines is in progress. 
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FIG. 1. 

7T-N scattering phase shifts. The dotted curve is the extrapolated 
HBChPT result, with the chiral parameters of (9) . The dashed line 
is the IAM with the same parameters, and the continuous line is 
the IAM constrained to resonance saturation (fit 2, see text). The 
dashed-dotted line is one unconstrained IAM fit. The data come 
from |nj. 

TABLE I. 

IAM results. In fit 1, we keep the 0(q 2 ) values of In fit 2, 
the a{ are constrained to the ranges predicted by resonance satu- 
ration j^] . Due to the strong correlations (see text) only certain bi 
combinations could be meaningful. 
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TABLE II. 

7rN threshold values with different IAM fits, and their experimental 
values (from Jlo| , see text). 
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